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Abstract 

We address some criticisms by Eletsky and loffe on the extension of QCD sum 
rules to finite temperature. We argue that this extension is possible, provided the 
Operator Product Expansion and QCD-hadron duality remain valid at non-zero 
temperature. We discuss evidence in support of this from QCD, and from the 
exactly solvable two- dimensional sigma model 0(N) in the large N limit, and the 
Schwinger model. 



Some time ago Bochkarev and Shaposnikov proposed an extension of the QCD sum 
rule program to non-zero temperature, and made an application to the two-point function 
involving the vector current. This application was reconsidered in 0. Later on we 
discussed the axial- vector channel |^ and the nucleon channel using Finite Energy QCD 
sum rules (FESR). The results from these analyses indicate a substantial rearrangement 
of the hadronic spectrum with increasing temperature, and hint on the existence of a 
deconfining phase transition. This was later confirmed in where a formalism valid 
even for T near the critical temperature was used. QCD sum rules Q are based on the 
Operator Product Expansion (OPE) of current correlators at short distances, suitably 
extended to include non-perturbative effects. The latter are parametrized in terms of a 
set of vacuum expectation values of the quark and gluon fields entering the basic QCD 
Lagrangian. Contact with the hadronic world of large distances is achieved by invoking 
the notion of QCD-hadron duality. The values of the vacuum condensates cannot be 
calculated analytically from first principles, as this would be tantamount to solving QCD 
exactly. Instead, they are extracted from certain channels where experimental information 
is available, e.g. e+e~ annihilation, and r decays 0. It is also possible, in principle, to 
estimate them numerically from lattice QCD. An extension of this formalism of QCD 
sum rules to finite temperature entails the assumptions that (a) the OPE continues to 
be valid, except that now the vacuum condensates will develop an (a-priori) unknown 
temperature dependence, and (b) the notion of QCD-hadron duality also remains valid. 
In analogy with the situation at T=0, the thermal behaviour of the vacuum condensates 
is not calculable analytically from first principles. Some model or approximation must 
be invoked, e.g. the dilute pion gas approximation, lattice QCD, etc.. The quark, the 
gluon, and the four-quark condensates at finite temperature have thus been estimated in 
the hterature [§- ig. 



In a recent paper, Eletsky and loffe [|TT] have criticized the QCD sum rule program 
proposed in and developed in In this note we wish to address this criticism, 

and hopefuly clarify the issue. We will argue that, provided the OPE and QCD-hadron 
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duality remain valid at finite temperature, the approach of |T]] is basically correct. 
We illustrate our argument with the current correlator involving the axial- vector current, 
although it can be trivially generalized to any local current operator. In addition, we 
provide further supportive evidence from two exactly solvable field theory models: the two- 
dimensional 0(N) sigma model in the large N limit, and the two-dimensional Schwinger 
model. 

The basic object to be considered is the retarded (advanced) two-point function after 
appropriate Gibbs averaging 

n (g, T) = z J d^x exp{iqx) 9{xo) « [Jix), j\0)] » , (1) 

where 

«A-B »= J2 exp(-E„/T) {n\A ■ B\n) /Tr{exp{-H/T)) , (2) 

n 

and I > is a complete set of eigenstates of the (QCD) Hamiltonian. The OPE of n(g, T) 
is formally written as 

n (g, T) = Ci « I » +J2 Cr{q) « Or » , (3) 

r 

where the Wilson coefficients Cr{q) depend on the Lorentz indices and quantum numbers 
of the external current J{x), and also of the local gauge-invariant operators Or built from 
the quark and gluon fields of QCD. The Cr{q) could also depend on the temperature, 
but since this is not essential for the argument we shall ignore this dependence. The 
unit operator / in Eq.(3) represents the purely perturbative piece. The OPE is assumed 
valid, even in the presence of non-perturbative effects, for < (spacelike), and | g^ |^ 
^QCD- principle, all Wilson coefficients are calculable in perturbative QCD to any 
desired order in the strong coupling constant. In the sequel we shall work at leading (one 
loop) order for simplicity. The non-perturbative effects are then buried in the vacuum 
condensates. Since these have dimensions, the associated Wilson coefficients fall off as 
inverse powers of = — g^. 
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For instance, if the current J{x) in Eq.(l) is identified with the axial- vector current 
='■ ^(a;)7^75rf(x) :, then with 

n^xi.(g,T) = -g^^Ui{q,T) + qf,q^Uo{q,T) , (4) 

one easily finds at T=0 [§] 

Stt no(g,T = 0) = -In— + + +■■■, (5) 

where /i is a renormalization scale, and e.g. the leading vacuum condensate is given by 

^4(04) = ^{as G^) - SA.iqq) , (6) 

with frig = {niu + md)/2, and < qq >=< uu >~< dd >. The function no(q'), Eq- (5), 
satisfies a dispersion relation 



defined in this case up to one subtraction constant, which can be disposed of by e.g. taking 
the first derivative with respect to in Eq.(7). The notion of QCD-hadron duality is 
implemented by calculating the left hand side of Eq.(7) in QCD through the OPE, and 
parametrizing the spectral function entering the right hand side in terms of hadronic 
resonances, followed by a hadronic continuum modelled by perturbative QCD. In this 
fashion one relates fundamental QCD parameters, such as quark masses, renormalization 
scales, vacuum condensates, etc., to hadronic parameters such as particle masses, widths, 
couplings, etc.. The convergence of the Hilbert transform, Eq.(7), may be improved by 
considering other integral kernels. This leads to other versions of QCD sum rules, such 
as the Laplace transform, FESR, etc.. 

After this introduction we address the criticism raised in [|ri|]. The states | n > entering 
Eq.(2) can be any complete set of states, e.g. hadronic states, quark-gluon basis, etc.. 
Eletsky and loffe claim that below the critical temperature Tc the suitable set of states 
is the hadronic set but not the quark-gluon basis. In support of this, they argue that 



for T < Tc the original particles forming the heat bath, being probed by the external 
currents, are hadrons. They go on to say that a summation over the quark-gluon basis in 
Eq.(2) would require taking into account the full range of their interaction, but that no 
account of confinement was given in []I| - 

First and foremost, we believe that the quark-gluon basis is indeed the appropriate basis 
to be used in QCD sum rule applications. If this were not the case, it would mean that the 
fruitful notion of duality would abruptly loose meaning as soon as the temperature is raised 
from T=0 to some arbitrary small value of T. This would be a rather bizarre scenario. 
According to the QCD sum rule philosophy, at T=0 one calculates the theoretical left hand 
side of Eq.(7) through the OPE, Eq.(3), i.e. one uses quark-gluon degrees of freedom, and 
duality relates this QCD part to a weighted average of the hadronic spectral function. 
The latter arises from using hadronic degrees of freedom. At very low temperatures 
the hadronic spectrum will change very little, and the external current will still convert 
into quark-antiquark pairs. Hence, it is only reasonable to assume that nothing drastic 
will happen with duality. At finite temperature, though, there is a new effect coming 
into play, i.e. there are contributions to the QCD and hadronic spectral functions in 
the space-like region (as opposed to only the time-like region at T=0). However, these 
additional contributions vanish smoothly as T approaches zero, i.e. they do not introduce 
any discontinuous behaviour. 

Besides, the fact that the heat bath is mainly composed of hadrons at small T is not 
in contradiction with the use of the quark- gluon basis. For instance, quarks enter the 
QCD perturbative term through loops, and have any value of momentum. One should 
emphasize here that the currents in (1) are external objects, and hence need not be in 
thermal equilibrium with the heat bath which is being probed. 



Next, contrary to the statement made in [11], confinement has indeed been taken into 



account in [|T| - [^] in the standard way, i.e. through the non-vanishing vacuum condensates 
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in the OPE, Eq.(3). The thermal beahaviour of these condensates is a separate matter. We 
have used in the chiral perturbation theory estimates of the thermal quark and gluon 



condensates |Q, and in |Q] we used the results of [g. Recently ||T0|, the T-dependence 
of four-quark condensates was obtained using soft pion techniques in conjunction with 
Eq.(2), where the summation was performed in the hadronic (pion) basis. In this instance 
we do agree that the hadronic basis is the appropriate one. However, this has nothing 
to do with QCD nor with duality . It is only one of many theoretical approximations to 
estimate the temperature dependence of the condensates. 



Finally, we wish to comment on an implicit claim of Eletsky and loffe |11|, which was 



stated more explicitly earlier in [|I2[ . This is, that calculating Gibbs averages in the quark 
basis at small T implies dealing with soft on-shell quarks, which are usually referred to 
as condensates. Although the quarks are on-shell, they can have any value of momentum 
when circulating around loops. And, there is no possibility of confusing perturbative 
contributions of on-shell quarks with quark condensates. This is obvious at T=0 if one 
computes the imaginary part of a current correlator, e.g. to one loop order, where the 
quark-antiquark intermediate state is on-mass-shell. For instance, the leading perturba- 
tive contribution to the imaginary part of no(g, T = 0) in the axial- vector channel, cal- 
culated using the Cutkosky rules (on-mass shell quarks understood) gives: Jmllo = l/Svr. 
This term cannot possibly be confused with the quark condensate contribution in Eq.(5) 
which, first, is real and, second, it has a different dependence. We argue next, that 
this is also the case at finite temperature. 

With = uj"^ — q^, and in the rest frame of the medium (q ^ 0), the imaginary part of 
no(ti;,T) in the axial-vector channel, to leading order in perturbative QCD, is 

i Im (c, q = 0, T) = ^ v{u^) [3 - v\u)]th (|,) e{uj' - Am]) , (8) 

in the time-like region, and 

-Imn(-)(u;,q = 0,T) = -i^5(a;2) T dzM^)[^ ' v\z)]2nJ^) , (9) 



in the space-like region, where v{x) = (1 —Arrig/x'^Y^'^ , and np is the Fermi thermal factor. 
The non-perturbative contributions to the OPE involve the quark and gluon condensates, 
both of dimension d=4, the four-quark condensate of dimension d=6, etc.. All these are 
real, and exhibit a temperature dependence very different from that in Eqs. (8)- (9). For 
instance, the quark condensate at low temperatures is of the form << qq >>=< qq > 
(1 — aT'^/f^), while the gluon condensate is essentially independent of T P|. In view of 
these differences, we see no possibility of confusing thermal quark loop contributions with 
thermal quark condensates. 

Additional evidence can be found in the framework of the 0(N) sigma model in the large N 
limit, and in the Schwinger model, both in two- dimensions. These models were used in the 
past |T^ in order to verify the validity of the OPE. Since these two models can be solved 
exactly, one can compare the results of the exact calculation with that from the OPE. They 
turn out, in fact, to be identical. We obtain below the thermal Green functions in these 
two models, and show that the temperature corrections to the perturbative contributions 
cannot be shifted to the non- perturbative terms, which develop their own T-dependence 
(calculable within the framework of these models). 

We consider first the 0(N) sigma model in 1-1-1 dimensions which is characterized by the 
Lagrangian 

C=^-[d,a%x)][d,a'^{x)], (10) 

where a = 1,...N and cr^cr'* = N/ f, with f being the coupling constant. In the large N limit 
this model can be solved exactly (for details see [|l^), it is known to be asymptotically 



free, i.e. 



= 1 t/ 2 (11) 

and in spite of the absence of mass parameters in Eq.(lO), it exhibits dynamical mass 
generation. In addition, in this model there are vacuum condensates, e.g. to leading 
order in 1/N 

(0|a|0) = v^m^ , (12) 
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whereas all other condensates factorize 



(0|a'=|0) = (VArm^)\ (13) 



In Eqs.(12)-(13) the a field is: a = f{d^a'^y' / yN . We are interested in the Green function 
associated with the propagation of quanta of the a field. At T=0, in Minkowski space, 
and at the one loop order, this Green function is given by 

(Pk i i 



r(Q) 



2 



(27r)2 P _ ^2 (g _ ^)2 _ ^ 

1 VQ^ - - 4m^ 

In , (14) 



^g2(g2 _ 4^2) VQ^ + VQ" - 

Expanding Eq.(14) at short distances leads to 

1 In(mVg^) / 2 (OHO) 6 (0|a^|0) \ 



A separate calculation, based on the OPE of the current correlator involving the scalar 
current Js = /((9^cr'*)^, gives exactly the same answer as Eq.(15) |p!3| . 



We have calculated the Green function Eq.(14) at finite temperature. Its imaginary part 
can be integrated analytically in closed form and is 

Imr(cu,q = 0,T) = -1^ [I + Sub (uj /2T)] 



1 

+ 2 



<< a » 6 << >> 



+ 



(16) 



where the first term above corresponds to the perturbative contribution, the second to 
the non-perturbative, and is the thermal Bose factor. Equation (16) is valid in the 
time-like region; the space-like region counterpart vanishes in 2 dimensions. Since the 
model is exactly solvable, the thermal behaviour of the vacuum condensates can also be 
calculated, viz. 

«a»=<a>[l + ^nB{uj/2T)] (17) 
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In this case the vacuum condensates contribute to the imaginary part, and as Eq.(16) 
shows, the thermal dependence of the perturbative piece cannot be absorbed into the 
condensates. Hence, no confusion should arise. We have not discussed above the effects 
of renormalization |T^, as they are not essential to our argument. 



Finally, we consider the Schwinger model in 1+1 dimensions, with the Lagrangian 



where = id^ + eA^. This model has been solved exactly |jT^, and in the framework 



of the OPE |T3]. The short distance expansion of the exact solution coincides with that 



from the OPE, as shown in |13|. Here, we are interested in the two-point functions 

H++(x) = (0|r{j+(x)j + (0)}|0) (19) 
H+_(x) = (0|r{j+(x)j-(0)}|0) (20) 

where the scalar currents are 

= i>R^L , r = i^R (21) 

with ipL,R = (1 ± 75)'^/2. The function H++(Q) vanishes identically in perturbation 
theory, and the leading non- perturbative contribution involves a four-fermion vacuum 
condensate. At T=0 one finds 



16 - 

H++(g,T = 0) = -g^(0|^ij d, ^R d, V-lIO) (22) 
On the other hand, the function H+_(Q) is purely perturbative. At T=0 it is given by 

m 



1 /, M2 \ 

n.-«.r = o) = - (i„^ + ^ + ...j , (23) 

where M is an ultraviolet cutoff. We have calculated the thermal behaviour of these 
current correlators and obtain, e.g. for their imaginary parts in the time-like region 
(there is no space-like contribution in 2 dimensions) 

ImH++(cu,q=0,T) =0 (24) 



Im n+_(cj, q = 0, T) = i [1 - 2n^(cu/2T)] (25) 

Hence, the choice of the fermion basis in the Gibbs average of current correlators does 
not imply confusing these fermions with condensates. The (perturbative) fermion loop 
terms and the (non-perturbative) vacuum condensates develop their own temperature 
dependence. In this particular example, non-perturbative terms are totally absent in 
/mn+_, and they do not contribute to /mn_|_+ because they are real. This peculiarity 
makes the argument particulary transparent. 

It should be mentioned, in closing, that there are some unresolved problems in finite 
temperature QCD sum rules, e.g. when one tries to use them to extract the thermal 



behaviour of the vacuum condensates ||T5[. However, these problems are unrelated to the 



arguments given by Eletsky and loffe [11 . 
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